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We study the properties of tlie trace T of tfie QED energy-momentum tensor in tfie presence of 
quasi-constant external electromagnetic fields. We exhibit the origin of T in the quantum nonlinear- 
ity of the electromagnetic theory. We obtain the quantum vacuum fluctuation-induced interaction 
of a particle with the field of a strongly magnetized compact stellar object. 



PACS numbers: 11.15.Tk,12.20.Ds,95.36.-|-x,97.60.Bw 



O 

(N 



Oh' 
CD 



> 
m 



00 
O 

> 

X 
J3 



Introduction 

Quantum electrodynamics (QED) in (quasi-)constant, 
homogeneous external electromagnetic (EM) fields pro- 
vides an opportunity to study the properties of the vac- 
uum state structure under the conditions of extreme 
external fields. In the presence of an electromagnetic 
field that varies negligibly on the space-time scale of the 
electron-positron fluctuations in the vacuum Ac = h/meC 
leads to an effective nonlinear electromagnetic th eory 
via the Euler-Heisenberg (EH) effective action |l|-[lO|. 
The physical observables and effective action induced by 
quasi-constant external electromagnetic fields are well- 
defined, because QED is an infrared-stable theory in 
which the electron mass me is the key scale parameter. 

The requirement of Lorentz symmetry admits only one 
essentially new contribution in the final expression for the 
energy-momentum tensor T^"^, a vacuum energy term 
proportional to g^'^ . This term is similar to Einstein's 
form of dark energy, the cosmological constant A. The 
related repulsive anti-gravity like effect of the energy of 
vacuum fluctuations may become accessible to labora- 
tory experiments: pulsed laser technology is advancing 
rapidly towards the 'critical' held strength 

Ec = Bc = mVe = 1.3 x 10^^ V/m = 4.4 x 10^ T, 

making strong-fleld QED effects possible in the labora- 
tory [11|. Magnetic fields of strength B ^ Be are also 
encountered in the stud y o f supernovae and post-main 
sequence stellar objects [ia - ll5l |. and we show here how 
the vacuum energy leads to a novel often repulsive force 
between the stellar magnetic field and a charged particle, 
influencing matter accretion and the supernova bounce. 

Much of what we present here is a general property 
of any nonlinear theory of electrodynamics [16|, appli- 
cable also to Born-Infeld elcctromagnctism, for exam- 
ple. We achieve considerable simpliflcation and insight 
exploiting a common feature of any nonlinear electro- 
magnetism, namely the presence of a dimensioned field 
scale Ec = m^/e which we express using a mass scale to. 
We write the (nonlinear) effective electromagnetic action 



Veff = -S + m /eff —7, 
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as a function of the (Lorentz) scalar S and pscudo scalar 
V. In Eq. (HI) /og contains solely contributions nonlin- 
ear in S, excluding a possible linear term, e.g. Slnm/fi 
where fi is another scale. When such term is included we 
omit the bar on fes. 



Energy-Momentum Trace 

We study the energy-momentum tensor for an effective 
nonlinear action (g — det g^ii^) 
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The dimensionless dielectric function e = —dVes/dS 1 
in the the classical Maxwell T^^x limit Kff —S- 
The form of Eq. © agrees with Eq. (A3) in [l^l and 
Eq. (4.17) in [l|]. 

The first term in Eq. ([S]) is traceless. The second term 
in Eq. ([5]) provides as noted the only possible covariant 
extension and is identically the trace of energy momen- 
tum tensor. Using Eq. ([T]) 
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Our separation of the trace T from the off-diagonal 
Maxwell-like part of the energy-momentum tensor iso- 
lates the field induced, gravitating energy-momentum of 
the vacuum. 

Before exploring the physical consequences, we must 
pause and clarify the precise meaning of T. Terms linear 
in the invariant S do not contribute to the right side of 
Eq. ([6]) since they cancel explicitly in the middle paren- 
theses. Only nonlinear (in S, V) EM-theories can have an 
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energy-momentum trace and the bar above VoS reminds 
us of this. This is another way to say that QED with 
massive electrons is ab initio not conformahy symmetric. 
Massive QED does not share the more challenging issues, 
such as conformal symmetry breaking, surroundi ng p ar- 
allel efforts in quantum chromodynamics (QCD) jl9l |. 

A better understanding of these remarks is achieved in 
QED by connecting T to the Diracf electron-positron) 
condensate induced in the vacuum [20| , which is directly 
related to the effective action 

— m(-0V') = mttr (Sp — S^) ~ m °^ . (7) 

dm 

The middle expression of Eq. ([7]) exhibits the condensate 
as the difference between normal ordering of operators 
in the no-field (also called perturbative) vacuum and the 
with-field vacuum. The argument of the derivative in 
Eq. ([7]) and in Eq. © is not exactly the same: The dif- 
ference in format between Eq. (jS]) and Eq. ([7]) is due to 
the leading linear term in S 
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Combining Eqs. (O, ©, and ([5]), we obtain the well 



known relation 12 ll 



T- 



2a 



2a2 



S + m^ipip) 



(9) 

4V' . J ' , (10) 

45toJ 

Eq. ([9|) displays two contributions to the trace of the 
energy-momentum tensor: gauge field and matter field 
fiuctuations. There is an important exact cancella- 
tion between these two terms. Were the first term in 
Eq. © (erroneously) omitted, the trace T would be that 
much greater, and for any applied magnetic field the 00- 
energy-density-component of the energy-momentum ten- 
sor would be negative. The QED vacuum would be un- 
stable, and the naive perturbative QED vacuum could 
reduce its energy by spontaneously generating a state 
with magnetic field. For all practical purposes, the form 
of Eq. © is confirmed by the observed stability of the 
QED vacuum and work claiming otherwise will need to 
address that important issue. 

In fact, the relative sign in Eq. ([9]) agrees with Eq. (5) of 
Ref. [l^ and Eqs.(35) in [l^l with the recognition that the 
(fermion) Gell-Man-Low /3-function in QED is positive 
definite. Our Eq. © (and the more explicit form of T, 
Eq. (|lip below) , agrees with [2^ , which result is a bit 
surprising since it follows from the clearly contradictory 
supposition that T = mdVcs/dm. Moreover, there are 
quite a few other instances in literature where the first 
term in Eq. ([9]) is omitted. 



Strong Fields Simulacrum of Dark Energy 

Applying Eq. ^ with the Eulcr-Hcisenbcrg-Schwingcr 
effective action (using Schwinger's notation and units in 



which a = e^/in 0) gives an explicit formula for the 
trace, 

G.e-'"^^ta6^^^-il (11) 



T = — <> 

Stt 47r2 7o 



tan(eas) 



wherein the invariant magnetic- and electric-like field 
strengths are 



reducing as indicated to the classical magnetic and elec- 
tric fields when one invariant vanishes. 

In numerical evaluation of the energy-momentum ten- 
sor for arbitrarily strong fields we employ the method 
developed in @. Consider first the stable field configu- 
ration B ^ Q,E = Q. The subtracted meromorphic (i.e. 
residue) expansions of the function 
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display the stabilizing change in sign following the second 
subtraction. The sums and integrals are absolutely con- 
vergent, so we may resum the resulting series, obtaining 
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in which /3' = -nm^/eB = -KjiBjE^). Eq. (US]) presents 
the (renormalized) effective action also seen in Ref. @. 
Numerical evaluations of the condensate Eq. ([M]) and 
trace Eq. ([15]) in the presence of magnetic and electrical 
field (real part only) are shown in figure [H We discuss 
elsewhere [i6[ the case that both E, B are non zero, how 
real and imaginary parts contribute together for electric 
fields, and the case of spin-0 matter fields. For fields way 
above critical the results presented in figure [T] are not 
accessible in practice by direct integration of the proper 
time representation Eq. PT|) 

For the electric field case, the corresponding meromor- 
phic expansions of x cot x show poles on the real s-axis. 
We assign to the mass a small imaginary component 
m? — 7> VP? -f ie, replacing in Eq. (|14p and Eq. ([T5|) the 
denominator z2 -f- 1 — > z2 _ ^ + te. Thus, in the presence 
of an electric field, there is also nonperturbative imagi- 
nary contribution to T ■ 

From now on in this work we address strong magnetic 
fields. T in the presence of a magnetic field is positive 
for any given field B, in contrast to the negative of the 
condensate m{tp'ip). The manifest signs of the two ex- 
pressions Eqs. p4p and psp . which determine the physics 
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FIG. 1: The condensate —m{tptp) (top) and the trace of 
energy-momentum tensor T^'' (bottom) in units of m'*, as 
a function of magnetic B (red) and electric E (blue, dashed 
and solid) field strengths. The negative of the electric field 
result is plotted where appropriate. The dotted (for E) and 
dashed-dotted (for B) lines show the weali-field expansions 
up to 



outcome of this investigation justify the time and effort 
spent showing how T does not include the term linear 
in S, while m{ip'ip) does. Clarification of this exclusion is 
necessary since as noted T and m{ip'il}) are often conflated 
in literature. 

The trace T gravitates, just as the Casimir energy 
does H^l- Because in the Euler-Heisenberg-Schwinger 
calculation the 'constant' external field is global in ex- 
tent, this energy- momentum is manifested in the form 
of a cosmological constant. In contrast to matter, for 
which the particle pressure acts outwards, the pressure 
part of energy-momentum tensor described by — T/4 acts 
inwards. This is a general feature of any 'false' vacuum 
state: the outside true vacuum the squeezes the false, 
higher energy density vacuum out of existence. 

The sign reversal of T pressure (compared to pressure 
of regular matter) overwhelms the gravity of the positive 
energy density, providing the anti-gravity effect associ- 
ated with dark energy. The similarity of T to the cos- 
mological constant was noted before by Schiitzhold [l^ 
and can be made explicit in the Einstein equation by 
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With a sign like that of dark energy K/AnG = A ~ 
(2.3 meV)'* = 4.1 x 10~'^^TOg, T is the dominant con- 
tribution in a domain of space with strong fields and is 
naively expected to generate a pressure that sweeps out 
matter, in analogy with the cosmological acting constant 
at large scales and pushing the universe apart. For com- 
parison, we note that the magnetic field-induced vacuum 
energy rivals cosmological dark energy A at Bd = 108 T, 
the scale of the largest static laboratory fields. 

Energy-momentum sources in Eq. (|16p involving local 
dark energy-like contributions have only recently been 
studied [24l - l27| . To see that the effects of local T are 
anti-gravitational just like cosmological A, we inspect the 
Oppcnheimer-Volkov equation 



dr c-^^ ° ^'r{r- 2GM) 



(17) 



where M{r) ~ f TQ^nr'^dr' and T^' = p — A is assumed 
isotropic. A does not contribute to the first term which 
is always non- negative: Tg -I- = {p + X) + {p — X), 
but A can make the second term {M + Airr^T^) change 
sign. Contributions to proportional to _g(^ (like A) thus 
weaken the pressure gradient and support heavier stars 
than otherwise expected. 

We have checked using Eq. ((T7)) that direct gravita- 
tional modifications to the mass-radius relation for com- 
pact stellar objects remain negligible, as might be ex- 
pected. At 60Bc, T = O.Sm^ = 1.14 x lO^^erg/cm^, 8 
orders of magnitude smaller than the pressures expected 
in the high density nuclear matter in the core of a post- 
main sequence star [l5| . However, this energy density 
is 2-4 orders of magnitude larger than the gravitational 
potential energy density of infalling stellar plasma. Thus 
while the insignificance of the QED T in gravity is a 
consequence of its objectively small energy density, the 
anti-gravitational effect of T suggests a closer study of 
the force experienced by individual particles in nonlinear 
electromagnetism is necessary. 



Particles in Overcritical Quasi-Constant Fields 

Forces present in the dynamical case can be much 
greater than those observed when the interacting bod- 
ies are studied in the hydrostatic equilibrium of Eq. (|17l) . 
The electromagnetic force determining individual par- 
ticle dynamics does not make relevant contributions, 
and the vacuum fluctuation-induced force contributes 
even less to Eq. (jl7p . Though very much smaller than 
Maxwell's force this force can be stronger than gravity 
and at times more relevant than the linear order force of 
Maxwellian electromagnetism. We will describe its fea- 
tures relevant to the charged particle dynamics within 
collapsing stellar objects. 
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Consider that the total ('t') electromagnetic energy- 
momentum tensor T^'^ due to both an external field 
('e') and a probe charge ('p') includes also an interac- 
tion energy- momentum Tj^jj , 



(18) 



with tensors T^'^ , T^'^ defined by the forms they take 
in isolation from each other. When the external field 
is much larger than the field of the probe particle, the 
electromagnetic energy-momentum tensor is expanded in 
the displacement tensor K^^'^ 



dV,s 



pfiv 



5/< 



cff 



(19) 



around the dominant contribution of the external field, 



p 



(20) 



with the subscript e reminding that derivatives are to be 
evaluated at the external field. 

The energy-momentum tensor T^"^ is expressed in 
Eq. ^ in terms of the field tensor F'^" , but only the 
displacement fields of the probe particle arc known ex- 
plicitly by solving Maxwell's equations with source 



By inverting Eq. (|19p we obtain 



dpt^^'dFap 



(21) 



(22) 



cff 



dFi^'^dF^s dF-r^dFa^ 



This rank 4 tensor transforms the functional dependence 
from the field tensor to the displacement tensor. We 
checked the validity of truncation by numerical evalu- 
ation of the derivatives of the action, which shows that 
(normalized) higher derivatives are suppressed even when 
the field is supercritical. A separable contribution of the 
Maxwell self-energy T^^" ~ ^Max p probe parti- 

cle is indeed found at next order d'^/dK'^ and would be 
subtracted. However, we find the effects of the terms 
from the second derivative are many orders of magni- 
tude smaller than those from the first derivative and do 
not discuss them further here. 
Using Eq. ^ in Eq. ^ gives 

Ttnt = T^p + (e - l)'T^c7 + a'^^^KpKf (23) 
where 

T^p = -(FrK^ + F::'K^n9.x + g'^ln^Kf (24) 

is the Maxwellian interaction with F^^K'^^ ~ 2(i?c • 

ifp - ^e'-Dp). The latter two terms of Eq. (^31) re- 
main after cancellation among the order a terms, and 



despite being order the last term (1 — e){dT/dS) = 
{d feff / dS){dT / OS) is kept for now. 

We view the net force (density) acting on the charged 
probe particle entering the domain of the external field 
in the usual way, requiring that inertial resistance bal- 
ance any breach of the conservation of the field energy- 
momentum, 



juFr + sr 



(25) 



using that —d,yT^p ~ j^F^" is the force obtained within 
Maxwell's linear electromagnetism. Here 



Sr~ie,~ ifjlF^- T^:d,{e,- if 



2-edT 



4 dS 



(26) 



with equality only approximate on account of the finite 
order expansions in T'"' Eq. ([20)) and a Eq. 



The conventional contributions of T^^p and F^^K^^ , 
i.e. the first term of Eq. ([23)) . to the force on the particle 
are obtained by integrating over a covariant hypersurface, 
and all frames being equivalent this integration is done 
most conveniently in the rest frame of the particle. 



d^x ElDi = - 



j (fx $c V'V^$p = S'^ J (fx $ePp 



(27) 

using that Dp = {Ze/r^)r (for a spherical charge) has 
only one component when choosing spherical coordinates. 

A physical charged particle has also a magnetic mo- 
ment fl, and thus a corresponding dipole magnetic field. 
As is well known this part of the force cannot usually 
compete with the effect of electrical particle charge. How- 
ever, in the present context we reach beyond the usual 
Lorentz force to the effect of vacuum fiuctuations and it 
is necessary to check if it is still justified to neglect the 
magnetic dipole in the strong magnetic field environment 
of a collapsing star. From the magnetic dipole interaction 
energy flp-Bc a force due to the gradient of the external 
magnetic field arises. Be changes on macroscopic scale 
though, and the parameter characterizing smallness of 
the effect is 1/mL ~ 10~^^ when varies on the scale 
L ~ 4 km. For comparison, the smallness parameter of 
the vacuum fluctuations arising from Euler-Hcisenbcrg 
action is {Be/ Bcfa/Abn. Seeing that vacuum fluctua- 
tion effects should dominate the magnetic dipole interac- 
tion for a stellar magnetic fleld B^ > 10~^ Be, we explore 
this domain further. 

Turning now to the latter two terms of Eq. ([23)) 
that represent the additional vacuum fluctuation-induced 
force, we observe that the gradient in the corresponding 
last term of Eq. ([25)) generates two contributions: the 
first as the gradient of Eq. ()27p and the second as the net 
change of the slowly-varying coefficient dT/dS over the 
domain of particle's field. The integrals over the parti- 
cle's field Eq. ()27p computed in the particle's rest frame. 
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the spatial components of the force Eq. f26|) on a point 
charge Pp = ZeS{x) in its rest frame are 



1 xr^ f< e-2ar 



E^e - Cc$eV5e, (28) 



where 



£-2a2r ise^ta ,,,, 



depends only on the scalar invariant of the external field 
Sc ■ On the right side of Eq. (|28p the gradient applied 
to the invariant S preserves the correct Lorcntz trans- 
formation property: although the potential $o appears, 
Eq. p8p is the gauge invariant correction to the linear 
force / = e{Ea + -y x B^) computed in the particle's rest 
frame. 

The component of Eq. (pS)) proportional to Co is qual- 
itatively different because, being proportional the gradi- 
ent ViS, it allows the transfer of energy from the magnetic 
field to in-falling particles. This property is in contrast 
to the first term in Eq. (|28p which produces a tiny change 
in the effective linear force (per mille at Be = -Be)- 

The weak-field expansion of Co Eq. (23) 




(30) 



obtained from the Euler-Heisenbcrg effective action, is 
usable up to ~ .1-Bc- The expansion shows that the 
predominant contribution to the gradient force is T, as 
the leading constant in Eg. (|30p is traced to d^T/dS^. 
Non-pcrturbative computation requiring employment of 
the nonperturbative numerical methods presented of the 
coefficients e — l,dT/dS,ctc. shows that the force per- 
sists in the considered high magnetic field domain despite 
what the perturbative expansion suggests. 



Particle Dynamics in a Stellar Magnetic Field 

As an application, we consider the dipole field of a 
strongly magnetized star. We are studying the force /r 
in the rest frame of the particle and in order to allow that 
the particle has a velocity relative to the star, we need 
to Lorentz-transform the field in the star's rest frame B' 
oriented at an angle V' from the direction of the parti- 
cle's motion. As a consequence of the transformation, 
the field of the star is seen by the particle to have an 
electric component. Specifically, i;^ = 7(1,0,0,/?) = 
(coshy, 0, 0, sinhy) and the Lorentz transformed fields 
are Bg = B'{cos ip z — j sin tp x) and Ec = —B'j(3 sin ip y. 
Using this in Eq. ([^5)) we obtain the force of the star's 
field on the moving particle. Note that (3, y can be posi- 
tive or negative. 

To compare with the gravitational force /g, we must 
also Lorentz transform it to the rest frame of the mov- 
ing particle. Although general relativistic corrections to 
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FIG. 2: The ratio /r//grav Eq. (fSTj) for a transverse {ip = 
7r/2) electron of rapidity y — cosh^^7 at the surface of a 
1.5M0, 12 km radius star with magnetic field given on the 
horizontal axis. Using the expansion Eq. (|30p for Cc produces 
the depicted values up to _B ~ .IBc. (Color online.) 



the Newtonian potential are significant near the stellar 
objects where such strong magnetic fields have been in- 
ferred, the Newtonian force is a reasonable first estimate 
modified only by multiplicative numerical factors of or- 
der unity down to a few times the radius of the future 
neutron star. 

We obtain the transformation property of the force /g 
considering the geodesic in the Schwarzschild metric (see 
e.g. Eq. 9.32 in (S^): transforming to the rest frame of a 
relativistic particle dilates the proper time, multiplying 
kinetic and total energies by 7^ . The energy equation for 
the geodesic therefore preserves its form if the same factor 
7^ is included also in the 'potential' terms, thus giving 
the transformation /g 1-^ The ratio of the radial 

vacuum fluctuation force to the Newtonian gravitational 
force for a transversely moving (ip = 7r/2) electron is 



/g 



3CcZel3jBsuii 



Rs 



eB. 



(31) 



illustrated in figure [H 

The ratio /r//g Eq. pip can be large as shown in fig- 
ure[5] due to the weakness of gravity, particularly for 
Be ^ 10~^Bc- The stellar magnetic field contribution to 
the gravitating energy density remains relatively small, 
yet it affects particle dynamics through the coupling to 
moving charge, through vacuum ffuctuation nonlincarity 
suppressed by lm\. Regarding sign in Eq. (|3ip . CPT 
symmetry of the vacuum assures that matter and anti- 
matter are expelled to the same degree: the effect of the 
force is the same for a left moving electron as for a right 
moving positron, as seen by simultaneously flipping the 
signs of Z and /3. Allowing for the distribution of charges 
and velocities with respect to the orientation of the field, 
we recognize that in a random medium (plasma) half of 
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the charged particles of each polarity at any given time 
is expelled. 

The interesting feature of the force Eq. ([28|) and 
Eq. pip is that a magnetic field gradient correlates 
velocity and charge. As the magnetic field curves tra- 
jectories of both left moving electrons and right moving 
positrons in the same direction, the noted symmetry 
in the effect of the vacuum fluctuation-induced force 
generates rotation in a net neutral plasma even while net 
current remains zero. Therefore, in a plasma made of 
negatively charged electrons and positively charged light 
nucleons, the matter which is ejected has opposite net 
momentum compared to the matter which is attracted 
and thus angular momentum is imparted to the magnetic 
source due to the mass asymmetry between positively 
and negatively charged particles. 

In summary, we have evaluated the trace T of the QED 
energy-momentum tensor and demonstrated that its gra- 
dient entails a significant and often repulsive force, which 
can be large compared to gravity, even while the rela- 
tive energy density of T remains small. Although the 
magnitude of the usual magnetic force v x Be is much 



larger than the vacuum-fluctuation induced correction 
Eq. (pS)) , only the latter is relevant in consideration of 
energy transfer to the particle and escape from the grav- 
itational potential well. The requisite energy exchange 
with a magnetic field, seen in the gradient ViS, is a conse- 
quence of the nonlinearity of the induced vacuum fluctu- 
ations, absent in classical Maxwellian electromagnetism. 

The quantitative study we present in Fig.[2]for the ratio 
Eq. (ini) indicates the force Eq. (gH) and Eq. (gS]), should 
have an impact on matter accretion and stellar collapse 
dynamics in astrophysical situations where strong mag- 
netic fields in excess of B ^ lO^^Bc = lO^T are known to 
exist. While treatment of the complete dynamical situa- 
tion is beyond the scope of this work, it is easy to imagine 
that this force could help the neutrino based transport 
phenomena [29| to propel the supernova bounce. 
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